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ASYMPTOTIC STABILITY OF DISSIPATED HAMILTON-POISSON 

SYSTEMS 
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Abstract 

We will further develop the study of the dissipation for a Hamilton-Poisson system introduced 
in [3] . We will give a tensorial form of this dissipation and show that it preserves the Hamiltonian 
function but not the Poisson geometry of the initial Hamilton-Poisson system. We will give pre- 
cise results about asymptotic stabilizability of the stable equilibria of the initial Hamilton-Poisson 
system. 
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^ 1 Introduction 

Q : Let 

^ ■ x = TlvH(x) 

^— » . 

be a Hamilton-Poisson system. If C € C°°(1R™,R) is a Casimir function then we have the following 

obvious equality II VC = 0. 

We will add to the Hamilton-Poisson system that we will denote by £n a dissipation term of the form 

GVC, where G is a certain symmetric matrix that will be discussed below and C is a Casimir function. 

The dissipated system £ will be 

x = ILvH(x) + GvC(x) (1.1) 

The notion of dissipative bracket was introduced by A. Kaufman [10] in his study of dissipative 
Hamilton-Poisson systems. In [15], P.J. Morrison introduced the notion of metriplectic systems which 
are Hamilton-Poisson systems with a dissipation of metric type. Metric type dissipation was also 
^— v . introduced in [15] and is given by a dissipative bracket constructed from a metric defined on the phase 

space. Dissipative terms and their implications for dynamics have been studied in connection with 
various dynamical systems derived from mathematical physics, see [8], [9], [11], [12], [13], [14]. 

The dissipation that we will study is the one introduced in [3]. Here we will write it in a tensorial 
form which enables us to better understand its geometry. It is a particular form of the metric dissi- 
pation introduced in [15]. More precisely, our system (1.1) is the same system with the one described 
by equations (25) in [15] under the condition (<? y )V-ff = 0. Like the dissipation from [7], our dissipa- 
tion conserves the energy. A dissipation form that conserves the symplectic leaves and dissipates the 
Hamiltonian function have been studied in [5]. 

We will also illustrate the improvement of the stability result obtained in [3] . This improvement is 
obtained by the observation that the dissipation conserves the energy. 

2 G = V#<g> V#- ||V#|| 2 I 

In this section we introduce the dissipation matrix 

G= S7H® VH -\\S7H\\ 2 I. 

This is the tensorial form of the dissipation constructed in [3] where the matrix I is the identity 
matrix on M. n . We denote by £ the dissipated system (1.1) in the case when the matrix G is as above. 



Lemma 2.1. We have the following properties for G: 
(i) G^H = 0; 

(ii) VC • GvC < 0, for any Casimir function C. 
The equality holds in a point x € R" iff VC(x) and V-ff (x) are linear dependent. 
Proof. For (i) we have 

GV# = (V#® V#)V#- ||Vi?|| 2 IViJ 

= (VF • VH)VH- \\VH\\ 2 VH = Q 

Analogous for (ii) we have 

VC-GVC = VC-[(Vi7®Vi7)VC-||ViJ|| 2 IVC] 
= VC-[(VC-VH)VH-\\VH\\ 2 VC} 

= (vc-vi7) 2 -||vc|| 2 -||vif|| 2 <o. 

by C-B-S inequality. Equality holds iff V-ff and VC are linear dependent. □ 

For the perturbed system £ the initial Hamiltonian remains a conservation low but the Casimir 
function C does not. 

Lemma 2.2. We have the following behavior: 

(i) H is conserved along the solutions oft;. 

(ii) C decreases along the solutions oft;. 

Proof. For (i) we have 

— h = x ■ vi/ = (nv# + GvC) ■ Vff 

= IIV# • VH + [{VH ® VH)vC - ||V#|| 2 VG] • VH 

since II is an antisymmetric matrix. 

We do the same type of computation for (ii) 

■4-c = x • vc = (nv# + gvc) • vc 

at 

= -vff-nvc + vc-Gvc 

= VC-GVC<0 
by Lemma 2.1 (ii). □ 

Next, we will study the relation between the set of equilibria E^ n for the unperturbed system £n 
and the set of equilibria E% for the perturbed system £. 

Proposition 2.3. We have: 

(i) If VC(i e ) 7^ for a point x e E W 1 then x e E E^ <^> VH(x e ) andVC{x e ) are linear dependent. 

(ii) E^dE iu . 



Proof. By definition we have that 

x e e E in O HvH(x e ) = 
and 

X e G E^ <=>■ nViJ(x e ) + GVC(> e ) = 

(i) Let x e G .Ef . If we multiply both sides of the above equality with VC(i e ) then 

VC(x e ) • nvi7(a; e ) + VC(x e ) • GvC{x e ) = 

«. -nvC(x e ) • vH{x e ) + VC(x e ) • GvC{x e ) = ^ VC(x e ) • GvC{x e ) = 

which by Lemma 2.1 (ii) implies that S/H(x e ) and VC(x e ) are linear dependent. 

For the converse, if VH(x e ) and VC(x e ) are linear dependent then there exists A G R such that 
VH(x e ) = AvC(x e ). Consequently, 

UvH(x e ) = AnvC(x e ) = 0. 

In the case VH(x e ) =0we obtain that G(x e ) = 0. 

In the case VH(x e ) ^ we obtain that VC(x e ) = jVH(x e ) and using Lemma 2.1 (i) we have 
Gs/C(x e ) = jGvH(x e ) = 0. In both cases we conclude that x e G E^. 

(ii) Let x e G E^. If VC(x e ) — then GvC(x e ) = and we obtain that HvH(x e ) — which implies 
x e e E^ u . 

If VC(x e ) ^ then there exists A G R such that vH(x e ) = AvC(x e ). We observe that LIV-ff(x e ) = 
AnvC(a; e ) = and consequently, x e G £^ n - □ 

3 Asymptotic stability 

We will briefly recall some definitions of stability for a dynamical system on R™ that will be used later 

x = f(x), (3.1) 

where / G C°°(R",R™). 

Definition 3.1. An equilibrium point x e is stable if for any small neighborhood U of x e there is a 
neighborhood V of x e , V <Z U such that if initially xq is in V , then (/>(t,Xo) G U for all t > 0. If in 
addition we have 



lim 4>{t,x ) = x e 

t— »oo 



then x e is called asymptotically stable. 



For studying more complicated asymptotic behavior we need to introduce the notion of w- limit set. 
Let </>* be the flow defined by equation (3.1). The w-limit set of x is uj(x) := {y G M. n \3ti,t2--- — > oo s.t. 
4>(tk, x) — > y as k — > oo}. The oj-limit sets have the following properties that we will use later. For more 
details, see [16]. 

(i) If <p(t, x) = y for some t G R, then ui(x) = ui(y). 
(ii) u)(x) is a closed subset and both positively and negatively invariant (contains complete orbits). 

The above properties of w-limit sets have been used in the proof of the following version of LaSalle 
theorem, see [3]. 



Theorem 3.2. Let x e be an equilibrium point of a dynamical system 

X = f(x) 

and U a small compact neighborhood of x e . Suppose there exists L : U — > R a C 1 differ entiable function 
with L(x) > for all x <E U \ {x e } 7 L{x e ) — and L{x) < 0. Let E be the set of all points in U where 
L{x) = 0. Let M be the largest invariant set in E. Then there exists a small neighborhood V of x e with 
V C U such that u(x) C M for all x eV. 

The next result describes the asymptotic behavior for the solutions of the dissipated system £. We 
introduce the following set C* = {x e M™ vC(x) = 0}. 

Theorem 3.3. Let x e s E% be an equilibrium point for the dissipated system £. Suppose there exists 
a function ip(H,C) E C°° (M 2 ,M) such that g^(H(x e ),C(x e )) > and x e is a strict relative minimum 
for L(x) = ip(H(x),C(x)) — ip(H(x e ),C(x e )). Then there exists a small compact neighborhood U of x e 
and another neighborhood V of x e with V C U such that every solution of £ starting in V approaches 
the largest invariant set M in U C\(E^ [J C*) as t — > oo (M is an attracting set). 



Proof. We have L(x e ) = with x e being a strict local minimum for L. Then there exists a small 

f)C [ 



compact neighborhood U of x e such that L(x) > and -A{H{x),C{x)) > on this neighborhood. We 



have the following computation 

= ^VH-UvH+^-VH-GvC+^C-mH+^vC-GvC 
on oH oC oC 

= ^-vh ■ nvff + ^-Gvh ■ vc - |^nvc . vii + % vc . Gvc 

an on oG oC 

Using Lemma 2.1 (i), antisymmetry of II and the fact that C is a Casimir function for II we obtain 
that 

L(x) = ^vC-GvC 
From the hypothesis tJj > on U and by Lemma 2.1 (ii) we have 

L(x) < 0. 

Using again Lemma 2.1 (ii) we obtain that E from the Theorem 3.2 equals E^ U C* which give us the 
result. □ 

Remark 3.4. Observe that L is also a Lyapunov function for the unperturbed system £n- Consequently, 
by adding the dissipation we render the stable points of £n into asymptotically stable points for £. Also 
the Casimir function C can be considered as a Lyapunov function. 

Corollary 3.5. Let x e e E^ and a function ip(H,C) € C°°(R 2 ,IR) such that §^(H(x e ),C(x e )) > 0. 
Suppose that the function L{x) = ip(H(x), C(x)) — tp(H(x e ), C(x e )) has the properties: 

i)SL{x e ) = 0; 

ii) S 2 L(x e ) is positive definite. 

Then there exists a small compact neighborhood U of x e and another neighborhood V of x e with 
V C U such that every solution o/£ starting in V approaches the largest invariant set M in U f](E^ [J C*) 
as t — y cxd. Moreover H remains constant along these solutions and C decreases along these solutions. 

Proof. It is easy to observe that the point x e is a strict relative minimum of L. All the condition of the 
Theorem 3.3 are satisfied and we obtain the desired result. □ 

Remark 3.6. If we consider, in the Corolarry 3.5, the function ip(H, C) = H + C we obtain the stability 
result of [3]. 



Corollary 3.7. In the hypotheses of the Theorem 3.3 we consider a point xo G V. Suppose that the set 
H~ l {{K\) P| U f](E^ (J C») has a unique point x^ (h = H(xq)). The solution x(t,xo) of the system £ 
which verifies the initial condition x(0,xo) = xq satisfies the property 

lini x{t 1 xo) = Xh- 

t— 7-OC 

We can interpret the above result as follows, using the Lyapunov function of Theorem 3.3 we obtain 
that every equilibrium point in the neighborhood V is asymptotically stable for the dynamical system 
on the corresponding level set. 

4 Applications to the rigid body dynamics 

The motion of a rigid body can be reduced to the translation of center of mass and rotation about it. 
Rotation is conveniently described, in a coordinate system with the origin at the center of mass and the 
axes along principal central axes of inertia, by Euler's equations. This equations can be written in the 
following form 

x\ = (j- - j-)x 2 x 3 + Ul 

X2 = (j- - -j^xix^ + «2 

X3 = (jj - 77)^1^2 +«3 

where x\ = J1W1, X2 = /2W2, %3 = Iz^z are the components of x, Ii,l2,l3 are the principal moments 
of inertia, u)i,u)2,u>3 are the components of the angular velocity and ui,U2,u 3 are the components of 
applied torques u. In this paper we suppose that I\ > I2 > I 3 . 

The system of free rotations, denoted by £n, has the property that u = 0. It has the following 
Hamilton-Poisson realization 

(( S o(3))*«R 3 , {-,.}_, H) 

where {• , •}_ is minus-Lie-Poisson structure on (so(3))* ~ M 3 generated by the matrix 

n = 




and the Hamiltonian H is given by 

I rf-^ rf^ rf^ 

H(x u x 2 ,X3) = ^ + | + |)- 

It is easy to see that the function Co G C°°(R 3 ,R) given by 

C (x 1 ,x 2 ,x 3 ) = -{x\ +xl +xj) 

is a Casimir of our configuration ((so(3))* w R 3 , {• , •}_), called the standard Casimir, i.e. 

{C ,f} = 

for each f G C°°(R 3 ,R) where 

{/,.9} = (V/) T n_V.g. 

The set of the Casimir functions is given by 

{C7 = ^(C )|^eC°°(R,R)}. 
It is known that the set of equilibria is given by 

E in = {(M , 0, 0) I M G R} |J{(0, M , 0) | M G R} |J{(0, 0, M ) | M G R} 



and the equilibrium points of the form (Mo, 0,0) or (0,0, Mo) are stable and the equilibrium points of 
the form (0, Mo, 0) with Mo ^ are not stable. 
In this paper we consider a torque of the form 

u = GvC, 



with C a Casimir function and G the matrix 

/ _ x 2 _ x 3 Xi'X 2 X\X 3 \ 

7J If hh 2 hh \ 

G = VH®VH-\\X?H\\ 2 1 = ^ -4-4 t 2 ^ 

11 " Jl-<2 It It J2J3 

2 2 

xix 3 x 2 x 3 _x t _ x 2 

\ hh hh If 7? 

For the case of the free rigid body the above matrix is the one also used in [15] 



Case I. We first consider the standard Casimir Co in order to construct the perturbation. We obtain 
a torque Uq with the components: 



*i[(* 



±\E2. 4_ (± 

h 1 1% ^ v h h> 1% 



h> h 



4 = *.[(£- £)% + (£ h^ 

%=^h-h) x T7+(h--hm 

Let £° be the rotation system with the torque uo . It is easy to see that the set of the equilibrium 
points of the dissipated system is 

E e =E in . 

Also we have 

C* -{(0,0,0)}. 

We consider Mo G M* . The point x e — (0, 0, Mo) is an equilibrium point of the rotation system £°. 
We define the function tp : M 2 — > R by 



M 2 



Cr, — I?,H. 



This function has the properties: 



and 



&(H.a,)-**-£) + i 



dtp 
dC~v 



(H(x e ),C {x e )) = l>0. 



Using our notations we introduce the Lyapunov function 



L(xi,X2,X 3 ) 



(C (xi,x 2 ,x 3 ) 



M 2 

-Y) 2 + C (x 1 ,x 2 ,x 3 ) 

1 M 2 1 

( f 2 , 2 , 2\ lvl \2 i 1 t 2 2 



I 3 H(x\,X2,X 3 ) 
x 2 ) _h [ d + f + f ) 

* H H ± 3 



It is easy to see that 



SL(x e ) = 



and 



S 2 L(x e ) = 



h 

7i 













1-f 

i2 






2M 2 



Using our hypotheses we observe that 5 2 L(x e ) is positive definite. Consequently, the hypotheses of the 
Corollary 3.5 are satisfied and we have the following stability result. 



Theorem 4.1. There exists a small compact neighborhood U of x e and an other V C U such that every 
solution of £° starting in V approaches U D {(0, 0, x) | X € K} as t — > oo. Moreover H remains constant 
along solutions and Co decreases along this solutions. 

If [/ is sufficiently small then U D -E^o C {(0,0, x) | sgn(x) = constant}. Suppose that sgn(x) = 1. 
Using the Corollary 3.7 we have the following asymptotic stability result. 

Theorem 4.2. I/xo eV i/ien tfte solution x(t,Xo) of the system £,o which verifies the initial condition 
x(0, xq) = xq satisfies the property 



lim x(t, x ) = (0, 0, v / 2I 3 H(x )). 

t— >oo 

We present the simulation of the rotation of the rigid body in the case when the principal moments 
of inertia are I\ = 4, Ii = 1.5, J3 = 1 and the initial conditions are (—0.1, 0.2, 0.175). 




Figure 1: The rotation of the rigid body 



Remark 4.3. We see that the function C\ = Co + (Co if-) 2 is also a Casimir function. If we 

consider ipi(H, C\) = C\ — I3H we have -^r- — 1 > and ip(H, Co) = tpi{H, C\). It is possible to apply 
Corollary 3.5 for the system £1 of the form (1.1) with the Casimir function C\ and the equilibrium point 
x e = (0,0,M o ). 



Case II. We study the case of the Casimir C2 = —Co We obtain a torque U2 with the components: 

1 = _a; i [(37 - 77)77 + (/7 ~ 77)77] 

Let £ 2 the rotation system with the torque 112. It is easy to see that the set of the equilibrium points 
of the dissipated system is 



E f 



E : 



in- 



Also we have 



C 2 * = {(0,0,0)}. 



We consider the equilibrium point x e = (Mo, 0, 0) with Mo £ 
If we define 



MH,C 2 ):=H+(C 2 
we have that 



2 ' + h ' 



*,„,,. M 2 l 
— (tf,C 2 )_2(C 2 + — ) + - 

and 

In this situation we can apply Corollary 3.5 and we obtain the stability result. 

Theorem 4.4. There exists a small compact neighborhood U of x e and an other V C U such that every 
solution of £ 2 starting in V approaches U D {(x,0,0) \x E R} as t -4 oo. Moreover H remains constant 
along solutions and C 2 = — Co decreases along this solutions. 

If C7 is sufficiently small then U fl £^2 C {(x, 0,0) | sgn{x) = constant}. Suppose that sgn{x) = 1. 
Using the Corollary 3.7 we obtain the asymptotic stability result. 

Theorem 4.5. If xq E V then the solution x(t,Xo) of the system £ 2 which verifies the initial condition 
x(0,a;o) = xo satisfies the property 

lim x(t,x ) = (y/2hH(x Q ),0,0). 
Case III. In [3] were introduced the functions 

c -ic M ° 2 i 2 Co -(r + M o\ 2 ,c 2 

Cs-(C -— ) - — _(C 2 + — ) + — 



and 



L 3 (a;) = F(x) + C 3 (») - fl"(» e ) - C 3 (x e ). 



It is easy to see that: 
i) SL 3 (x e ) = 0; 
ii) <5 2 i3(x e ) is positive definite because the Hessian matrix of L 3 in x e is 



M 2 











1 1 











1 



We denote by £3 the system of the form (1.1) with the Casimir C3. In [3] has been proved the 
following result. 

Theorem 4.6. There exists a small compact neighborhood U of x e and an other V d U such that every 
solution of £ 3 starting in V approaches U D {(x, 0, 0) | x £ M.} as t — > 00. 

Our results applies to equilibria that are already stable. By adding the dissipation we can make 
them asymptotically stable if certain conditions are satisfied. Going from asymptotic stability on the 
level set to stability in the whole space have been studied in [2] . By adding other type of controls one 
can stabilize the unstable equilibria. This subject is studied for example in [1], [4], [6], [17]. 
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